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We address at the mean field level the emergence of a Pomeranchuk instability in a uniform 
Fermi liquid with central particle-particle interactions. We find that Pomeranchuk instabilities 
with all symmetries except I = 1 can take place if the interaction is repulsive and has a finite 
range ro of the order of the inter-particle distance. We demonstrate this by solving the mean field 
equations analytically for an explicit model interaction, as well as numerical results for more realistic 
potentials. We find in addition to the Pomeranchuk instability other, subtler phase transitions in 
which the Fermi surface changes topology without rotational symmetry-breaking. We argue that 
such interaction-driven topological transitions may be as generic to such systems as the Pomeranchuk 
instability. 



I. INTRODUCTION 

Experimental evidence of "hidden" phases of itin- 
erant electron systenisi*2i^ and the prospect of real- 
izing novel conditions in layered heterostructures and 
ultra-cold gases have led to increased efforts to iden- 
tify unconventional phase transitions and predict their 
manifestations. To give three examples: the Fulde- 
Ferrell-Larkin-Ovchinnikov state has been proposed in 
organic superconductora^, superconductor-ferromagnet 
heterostructures^ and in imbalanced mixtures of ul- 
tracold atomsiSii; a supersolid phase is a possibility 
in Bose gases loaded on optical lattices^; and a "d- 
density wave" may be realized in ladder compoundsiSiii 
and possibly "hide" in the phase diagram of cuprate 
superconductorsSi, where other hidden order parameters 
have been proposed^2i^. 

In this context there has been a surge of interest in 
the Pomeranchuk Instabihty (PI)i2i. Through it a Fermi 
liquid may enter a "nematic" state characterized by a de- 
formed Fermi surface. It has been argued that such an 
instability may take place in quantum Hall systemsiii^ 
and in the metamagnets Sr3Ru2072, and URu2Si2i^. 
Moreover there is evidence that the Hubbard model has 
a phase with a distorted Fermi surfaceiSiiLiSiiSiSSiSi and 
the Emery model of a Cu02 plane has been shown to have 
a nematic ground state in the strong coupling limiljSS. 

More generally the PI is an interesting candidate un- 
conventional phase transition on account of its sub- 
tlety. Thus, considerable effort is going into charac- 
terizing it theoretically^'^i^'^i'^^i^'^i^^i^^i'^^i^" on the basis 
of phenomenological models featuring anisotropic effec- 
tive interactions. This approach is proving very suc- 
cessful in establishing some generic features of the phase 
diagram2i2iiS& and describing collective excitations and 
quantum critical fluctuation92SiSL22i22i2£. On the other 
hand it sidelines the question of how the anisotropy 
emerges in the first place^^ and what other^^-, perhaps 
even subtler instabilities may generically arise in such 
contexts. It is these questions that we address here. 



In this paper we present a mean field (MF) theory 
of the PI in a three-dimensional, uniform fcrmion liquid 
with a central effective interaction potential V{r). The 
authors of Ref. have pointed out that such interac- 
tion may lead to a PI. Here we show that the emergence 
of the anisotropic state from a Galilean invariant fluid 
requires repulsion with an intermediate range of the or- 
der of the inter-particle distance. This is confirmed by 
explicit calculation for a model interaction potential for 
which the theory can be solved analytically. However 
we also find that the intermediate-range repulsion leads, 
quite generally, to a different instability in which there 
is no symmetry breaking but the topology of the Fermi 
surface changes. We discuss the nature of this subtler 
quantum phase transition. A few instances of the two 
distinct types of Fermi surface shape instabilities that we 
find are pictured in Fig. ^ These two types of instability 
compete and we show that this conclusion is robust when 
we consider more more realistic finite range interactions. 



II. MEAN FIELD THEORY 

To motivate a microscopic theory of the PI we start by 
recalling the original, phenomenological theory due to 
Pomeranchuk^^;. Like him, we start with an unpolarized 
Fermi sphere and consider an infinitesimal change of the 
occupation numbers, N\^^„„ — >• A^k, o-o- + 5N\^^„„, arising 
from an angle-dependent modulation of the Fermi vector, 
kp — > kp + Skp (0), kp — > kp + aSkp (0) in the symmetric 
or antisymmetric spin channel, respectively^^. We then 
use Landau's expression for the corresponding change in 
the ground state energy: E ^ E + 6E, with 

6E = ^£(k)(5^k (1) 
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Figure 1: (color online) Different shapes and topologies of 
the Fermi surface. From left to right, in the first row; (a) 
unpolarized, undeformed Fermi sphere; (b) Fermi surface with 
an I = 2 Pomeranchuk deformation; (c) 1 — 3. In the second 
row; (d) Fermi sphere surrounded by an additional sheet of 
occupied states; (e) with a "hole" of vacated states at the 
center; and (f) with a shell of vacated states. 



Here SN^, = EaSNk.aa and SSI = 5 E.,^ '^^^'^^k,^.- 
Requiring 5E < leads to the PI conditions 



l + i^°'7(2? + l) <0, 
in terms of the Landau parameters, defined hy^ 
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where fc^^ is the radius of the Fermi sphere, Ui? is the 
Fermi velocity and Pi (x) is the Z"^ Legendre polynomial. 
For I = 0, Eq. (jSJl describes a quantum gas-liquid tran- 
sition (in the symmetric spin channel, s) or a Stoner in- 
stability (in the antisymmetric channel, a). For Z > it 
describes a Pomeranchuk instability. 

One crucial aspect of Pomeranchuk's theory is that it 
describes the instability in terms of the phenomenological 
Landau parameters, Fi'"". Here we want to establish the 
mechanism whereby the PI could take place in a system 
with a given microscopic Hamiltonian of the form 



(4) 
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where y(|r — r'|) is a local, non-retarded, spin- 
independent and central interaction potential. We ad- 
dress this question using MF theory. Our MF Hamilto- 
nian is 



Ho = (k)c|,_^Ck,^, 



(5) 



where cj^^ = fi"^/^/ c?-°re^*'' ''c^g.. This describes inde- 
pendent electrons with an arbitrary dispersion relation 
Ea (k) , which we treat as our variational parameter. Note 
that this MF couples only to the occupation number in 
k-space, Ak.o-o- = C|{j.^Ck,CT. Our theory thus preserves 
translational and gauge symmetry, but it can neverthe- 
less break rotational symmetry if the dispersion relation 
becomes anisotropic. For example, a nematic Fermi liq- 
uid state may be entered through a PI. It is an example 
of an "electronic liquid crystal state"—. 

Although our main results refer to the ground state, 
the derivation of the basic equations of the theory is 
much simpler at finite temperature. We thus approx- 
imate the free energy by F w {H — Ho)q + Fq, where 
(. . .)o = Z^^Tr {e-^^o . . .} with Zq = Tr {e-'^^o} ^nd 
Fq — —(3^^ \\iZq. It takes the form 



F = ^|7Vk,4-;^E^(|k-k'l)iVk^. (6) 
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where V = J (PHV (|R|) is the uniform component of the 
interaction potential and V (K) = J tfRe~^'^-'^V {\R\) 
its Fourier transform. The occupation numbers in k- 
space are given by 



(7) 



Requiring that F be stationary yields 

(k) = ^^-^,+ 1 S,,,,Vi\k - k'l)} iVk.,... 
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In the low-temperature limit, /3 — > oo, Eqs. (|6I7|I become 
E = J2N^,A-l-Y,V{\k-k'\)N^,,, (9) 
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and our MF theory is equivalent to trying variationally 
the following ground state: 
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Eq. jnj is our MF approximation to the ground state 
energy. Variation with respect to the occupation num- 
bers yields an expression identical to Eq. J^l, except that 
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the phenomenological functions Sa (k) (which here may 
depend on the spin) and /*'° (k, k') are now derived from 
our microscopic parameters via Eq. and 



r''^(k,k') = ^ jryy-ii^(lk-k'i) 



(12) 



where 7y = 1, in the s, a channels, respectively. 

Let us pause briefly to note the following subtlety. 
Eqs. (|9I10|I . from which all the subsequent results fol- 
low, could have been derived by minimizing the energy of 
the trial state given in Eq. Ijlip . However note this only 
determines the Fermi surface, but it under-determines 
the dispersion relation ecr(k). The justification of the 
particular form given in Eq. Q thus relies on the as- 
sumption that the low-lying excitations correspond to re- 
arrangements of the electrons in momentum space, whose 
energy is given by Eq. ^ [or, equivalently, that the equi- 
librium state at finite temperatures can be adequately de- 
scribed by the mean field Hamiltonian of Eq. Q .] This 
is necessary to justify the language we use below, e.g. in 
defining the Fermi velocity in terms of ecr(k). We stress, 
however, that the results themselves refer only to the 
equilibrium shape of the Fermi surface in the ! ground 
state and are therefore more general, and independent of 
the meaning assigned to ecr(k). 

To study the PI in our microscopic model we postulate 
an unpolarized, spherical Fermi surface 

Nk^a ^eikp- |k|) , (13) 

completely described by the Fermi vector kp > [Fig.Q] 
(a)], and use the above equations to determine whether 
the system has a PI. In the state described by Eq. (|13|l . 
the electron dispersion relation of Eq. © is given by 

/ \ 91-2 /-oo 

s{\k]) = |-(|kr-4)-^y^ drrVir)iU) 

ji {kpr) [jo (|k| r) - jo (/ci^r)] . 

This yields the following expression for the Fermi veloc- 
ity: 

VP = -kp + —£^ / dr V (r) ji (kprf (15) 
m hn Jo 
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where in the second line we have expressed vp m terms of 
one of the couping constants defined by Eq. ((TH|l . below. 
Note that the state described by Eq. ifT^ requires vp > Q. 
Together Eqs. (|12|l and (|15|l give the Landau parameters 
in Eq. ®, 



(27r)3 hvp 



(16) 



in terms of the microscopic parameters of the model. 
This, in turn, allows us to express the PI equations as 



Vi - 7^6i.oSt:V > ^^^vp, 



i2Trhf 
mpp 



(17) 



where the strength of the interaction potential in a given 
angular momentum channel Z = 0, 1, 2, ... is given by 



Vi = (Any / drr'^V (r) ji {kprf 



(18) 



For antisymmetric instabilities, 77 = 0, Eq. H17() takes the 
following, explicit form: 



dr Attt^V [r) ji{kpr) — ji{kpry 



> 



mkp 



(19) 

Eq. H17|l is our microscopic expression of the PI con- 
dition of Eq. (O. It is valid, within our MF ansatz 
of Eq. (O, for any system whose Hamiltonian has the 
form given by Eq. From it we can derive a series of 
conclusions concerning a Pomeranchuk instability in an 
isotropic system with central interactions: 

1. Purely attractive interactions can only lead to the 
gas-liquid transition {I — 0,7] — ly^; conversely, 
purely repulsive interactions can only lead to the 
Stoner or PI. 

2. There are no / = 1 Pomeranchuk instabilities. This 
is the type of pii^j^LS^ where rotational symmetry- 
breaking is achieved by displacing the Fermi surface 
so as to set up a charge (s) or spin (a) current. This 
is quite a general consequence of the well known 
relation between the effective mass and the Lan- 
dau parameter, Ff, in a Galilean invariant system, 
which is captured by Eqs. (|15|l and H16|l . On the 
other hand for spin-dependent interactions (or in 
lattice systems) , not considered here, we may have 
Ff 7^ Ff and then the instabilities considered in 
j^gfgj^jSijSe could be reahzed. 

3. The PI for / > 2 is degenerate in the spin chan- 
nel. These instabilities break rotational symmetry 
by changing the shape of the Fermi surface, with- 
out generating any currents of charge or spin — see 
Fig. njb,c). Our result implies that, at the insta- 
bility, it does not matter whether the lobes of the 
spin-up and spin-down Fermi surface point in the 
same direction^. Note this is quite different from 
the situation at I = (see point ^ above). 

4. Finally, from Eq. (|18|l we can also deduce that 
Vi — Vi cannot be large and positive, as required by 
Eq. (|17|l . if the repulsive part of the interaction is of 
very short range ro <?C kp^. In effect, r^ji {kpr)^ ~ 
,-2(i-i-i) Jqj. j, ^ ^-1 gQ Jqj. gm;]-^ short-ranged inter- 
actions Eq. lUHJ) gives Vi - " drr^^+i) - rg'+^ 
whence for small ro Eq. (|17|l can only be satisfied 
for / = 0. The extreme case of this is the repulsive 
contact potential V^(|r|) = |u| (J^'^) (r) , for which 
ro — Q and Vi = At: \u\ Si^. For this potential, our 
theory leads only to the Stoner instability. 
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Figure 2: (color online) Phase diagram showing zero- 
temperature instabilities of the types depicted in Fig. for 
the repulsive delta-shell model. The Stoner instability oc- 
curs above the solid line. The other instabilities described 
in Fig. occur above the long dash line (b), medium dash 
line (c), short dash line (d), long dash-dot line (e) and short 
dash-dot line (f). The colors correspond with Fig0 



III. DELTA-SHELL MODEL 

From point 01 we conclude that repulsive interactions 
with range at least of the order of the Fermi wavelength, 
ro > kp^, are necessary for the PI in an isotropic system 
with isotropic interactions. We investigate this further 
by choosing a specific form of the central interaction po- 
tential, namely the "delta-shell" potential: 

F(|r|)=g<5W(|r|-ro). (20) 

This is an idealization of an interaction with a very sharp 
peak at a particular inter-particle distance, |r| — tq. The 
"coupling constant" g has dimensions of energy x length 
and represents the product of the height and width of the 
potential barrier. 

For g < 0, this interaction potential can lead to super- 
conductivity with unconventional pairing22i2£. Likewise, 
we expect that for 5 > it will lead to a PI. Indeed 
Eq. gives 

Vi^iAirfgr^jiikprof (21) 

so, depending on the value of kpro, any value of / may 
become dominant. 

The particularly simple form of the interaction poten- 
tial in Eq. (|20|l allows us to write the key expressions in 
our theory of the PI analytically. In particular Eq. I|15() 
giving the Fermi velocity on the Fermi sphere reads 

h 2{kprnf . , .2 

vp = —kp + g ji [kprQ) . (22) 

m fm 

Together with Eq. H21|l and V = g^i^r^ these equations 
allow us to write the following, simple form of the Stoner 



and PI equations: 

ji {kprof - ji {kprof > j ' ^ = 0, 1, 2, ... 

zmrQ gkp 

(23) 

If the interaction is very strong, gkp ^ H'^/2mrQ, this 
gives a sequence of fixed phase boundaries at ji (kpro) = 
iji {kpfo) ■ In the opposite limit of very weak interaction 
the unpolarized Fermi sphere is, as expected, stable. 

The solid, long dash and medium dash lines of Fig. [21 
are the phase diagram obtained by solving Eq. (|^ for 
I < 3. Note that there are effectively only two, di- 
mensionless parameters in the theory22i2£: the "effec- 
tive range" kpro and "coupling constant" g/roeo (where 
So = /2m.r^). The first of these parameters is the range 
of the interaction measured in units of l/kp. The sec- 
ond is the product of the width of the potential barrier 
measured in units of that range and its height measured 
in units of the corresponding "localization energy" Eq- 
As expected, for small range, kpro < 3, we only find 
the Stoner instability. For longer ranges or, equivalently, 
higher densities, and sufficiently large values of the di- 
mensionless coupling constant (which we note that de- 
pends not only on g but also on m and ro), the PI can 
take place. 

IV. TOPOLOGICAL TRANSITIONS 

In addition to the anticipated PI, our analytic treat- 
ment of the delta-shell model also reveals a competing 
class of Fermi surface instabilities which, unlike the PI, 
occur without symmetry breaking. Consider the electron 
dispersion relation in the isotropic state, Eq. 114|) . For 
the delta-shell potential it takes the form 

£ (|k|) = 1^ (|k|^ - 4) - (kpro) (24) 

bo (|k|ro) - jo (kpro)] . 

This is plotted in Fig.Ofor three different values of kpro- 
The free-electron dispersion relation is modified by an os- 
cillatory term due to electron-electron interactions. The 
period of the oscillations is ~ r^^ . For small g/roEo the 
effect of these is the usual renormalization of the effec- 
tive mass m* = pp/vp, which follows from Eq. H22|l . 
However at large g/roSo the effect of interaction on this 
"bare" dispersion relation cannot be described simply as 
a renormalization of m. In fact it can lead to a dramatic 
change of the state of the system as the amplitude of the 
oscillations becomes large enough that either 

1. the dispersion relation dips below the Fermi level 
somewhere outside the Fermi sphere [Fig. 121(a)], 

2. it goes above the Fermi level at the center of the 
Fermi sphere [Fig. (31 (b)] or 

3. it peaks above the Fermi level at some intermediate 
k,0<k<kp [Fig. 131(c)]. 
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Figure 3: Electron dispersion relation for the 
"delta-shell" model in the state with an unpo- 
larized Fermi sphere: (a) kpro = 4, g/roso — 
(solid), 20 (long dash), 97.66 (short dash), 150 (dotted); 
(b) kpro — 6, g/roEo = 0, 2, 8.84, 16 (same order); (c) 
kpro = 9, g/roeo = 0, 20, 41.40, 60. 



In either case, Eq. (|10|) no longer reduces to Eq. (|13|l . 
Instead, either a thin shell of occupied states forms out- 
side the Fermi sphere [Fig. ^ (d)] , or states inside the 
Fermi sphere become vacated [Fig. ^ (e,f)]. The associ- 
ated instabilities are quite distinct from the Stoner and 
PI, as the change JA^k.o-, although infinitesimal, takes 
place away from the Fermi surface. Instead, they are 
continuous phase transitions in which no symmetry is 
broken, but the topology of the Fermi surface changes. 
In that sense they are more reminiscent of the Lifshitz 




Figure 4: Evolution of the Fermi vector as a function of the 
dimensionless coupling constant g/roeo for fixed particle den- 
sity N/fl — 9 r^'^. The plot shows both the "original" Fermi 
vector kp and the "emerging" Fermi vector k'p, namely the 
radius of the sphere of empty states depicted in Fig. 0(e), as 
the system enters from below the corresponding dome in the 
phase diagram (see Fig.l^J. 



transition22il£. There is, however, a crucial difference, 
namely that the present phase transitions are driven by 
electron-electron interactions, which induce the fermions 
to "migrate" to other regions of reciprocal space, rather 
than by the band structure. On the basis of this one 
would expect the present instabilities to have a much 
stronger thermodynamic signature. For example, the 
transitions illustrated in Fig.^d) and (f), with their un- 
derlying dispersions of Fig. |3Ia) and (c) respectively, re- 
sult in the appearance of entire new Fermi surface sheets 
with finite kp- This will lead in mean- field theory to a 
discontinuous jump in the density of states and hence in 
C/T. 

A general framework to understand such topologi- 
cal quantum phase transitions has been put forward in 
Ref. 113. In this formalism, the Fermi surface is a vor- 
tex "loop" in a four-dimensional space and the phase 
transitions we have just described correspond to the nu- 
cleation of new loops. One can also view these insta- 
bilities as generalizations to dimension larger than one 
of the phenomenon of "quantum Hall edge reconstruc- 
tion". The latter can be described as the emergence, 
due to interactions, of new Fermi points in the one- 
dimensional chiral Fermi liquid on the edge of a quan- 
tum Hall system^242ii4iiS-. Indeed in the instabilities de- 
scribed here always one of the new Fermi surfaces has 
negative Fermi velocity - analogous to the creation of 
left-moving quasiparticles in a right-moving chiral Fermi 
liquid™. Yang and Sachdev'*^ have recently described 
the quantum critical fluctuations for a phase transition 
of type HI above [Figs.[T](d) and El (a)]. 

It is important to note that the plots in Fig. O corre- 
spond to evaluating Eq. (|24|l at a fixed value oi kp- For 
a fixed number of particles, N, such solutions are valid 
only up to the instability, as beyond it they would violate 
Luttinger's theorem. To describe the migration of elec- 



6 




Figure 5: (color online) Phase diagram showing zero- 
temperature instabilities of the unpolarized, spherical Fermi 
surface for the "hard-core" model. The coupling constant is 
measured in units of the Fermi energy ep = h^kp/2m. Key 
as in Fig|5| 



trons in reciprocal space mentioned above, which hap- 
pens beyond the instability, it is necessary to determine 
kp self-consistently by requiring that the total number of 
particles be fixed. For example, for an instability of the 
type 121 above [Figs.|31(b) andn](e)], this means that 



N 



1 An 



(25) 



where k'p is the "emerging" Fermi vector at the centerer 
of the Fermi sphere, determined by e{k'p) — 0. This is 
demonstrated in Fig. 01 Note that the way kp grows as 
a function of coupling suggests thinking of this quantity 
as a sort of topological "order parameter" . 

It is evident from Fig.|51that the new sheet of electron- 
like or hole-like Fermi surfaces are initially formed by 
localized states, with Wi? = 0. As we progress into the 
new state, vp becomes finite. Conversely, if we run the 
process backward, the effective masses on the additional 
Fermi surfaces diverge (except for the hole-like Fermi sur- 
face in Fig.n](e), for which the Fermi vector goes to zero 
at the same time as the Fermi velocity). 

The short dash, long dash-dot and short dash-dot lines 
on Fig. 121 show the boundaries of these "Fermi surface 
topology transitions". Notably, except for a range of 
densities near where the / = 1 PI would have been, the 
unpolarized Fermi sphere is only stable at small coupling 



V. HARD-CORE MODEL 

Our results suggest that the essential ingredient of 
the PI in the present context, namely the finite range 
''0 ~ kp^, also leads to the interaction-driven Lifshitz 
transition. To probe the generality of this observation, 
we have repeated the calculation (this time by evaluating 



the mean-field equations H14() and (|17|l numerically) for 
a repulsive "hard core" potential 



V{\r\)^VQ (ro-|r|). 



(26) 



We have again found that, for tq > fc^^, there are, in 
addition to the Stoner instability, PI with / = 2, 4, 6, . . . 
Moreover we also find an interaction-driven Lifshitz tran- 
sition which, for certain ranges of values of kpro, takes 
place before the Stoner or PI set in. A phase diagram is 
presented in Fig. 

Unlike the delta-shell potential, for the hard core po- 
tential the PI domes are contained within the Stoner 
ones, i.e. the PI can only take place on a polarized Fermi 
surface, or for spinless fermions. However note that there 
are regions of the phase diagram where the boundary of 
the / — 2 instability nearly overlaps with that for Z = 0, 
indicating that the two transitions happen almost simul- 
taneously. 

The plot does not show all the instabilities: the I > 
2 PI take place at higher values of kpro than those 
shown. The are also other domes of topological insta- 
bility, though for this potential all of them are of the 
type in Fig.l (e). 

Our results for the hard-core potential not only sup- 
port our identification of a sharp feature of V (r) at 
fo^kp^ as the crucial ingredient for a PI, but suggest as 
well that, in such situations, the interaction-driven Lif- 
shitz transition is at least as likely to occur as the Stoner 
or PI. 

We conclude this section by noting that a similar anal- 
ysis using the screened Coulomb interaction 



V{\v\) 



2 

e e 



r|/ro 



(27) 



does not reveal either PI or instabilities of the Fermi sur- 
face topology — only the Stoner instability is realized, 
and that only if we allow rg to deviate from its Thomas- 
Fermi value. 



VI. CONCLUSIONS 

In summary we have studied, at the mean field level, 
PI of a uniform Fermi liquid with central fermion-fermion 
interactions. We find that PI of different symmetries may 
emerge from repulsive interactions of sufficiently long, but 
finite range tq ^ kp^ (with the interesting exception of 
the Z = 1 PI which never takes place). We have confirmed 
this by solving the theory analytically for an explicit form 
of the interaction potential featuring repulsion at a par- 
ticular distance tq. Surprisingly we have found that, in 
addition to the PI, there is also a new type of Fermi 
surface instability: the interaction-driven Lifshitz tran- 
sition. This topological phase transition is even subtler 
than the PI and seems to be generically associated with 
the class of models leading to the PI. Further support 
for this picture is provided by analysis of an additional 
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model, featuring hard-core repulsion. On the other hand 
the screened Coulomb interaction does not lead to these 
effects suggesting that a sharp feature (either a spike of 
repulsion or a sudden drop) must be present at the dis- 
tance Vq. 

Unlike the Lifshitz transition, the new quantum phase 
transition that we have described is fundamentally driven 
by interactions. Thus one would expect it to have a 
stronger thermodynamic signature. It will be of great 
interest, in the near future, to establish this signature 
and the properties of the novel state of matter this phase 
transition may lead to. 
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